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CN Abstract 

P-i| We prove the global existence of weak solutions to the Navier-Stokes equations of compressible 
" heat-conducting fluids in two spatial dimensions with initial data and external forces which are 
large and spherically symmetric. The solutions will be obtained as the limit of the approximate 
solutions in an annular domain. We first derive a number of regularity results on the approximate 
physical quantities in the "fluid region" , as well as the new uniform integrability of the velocity 
and temperature in the entire space-time domain by exploiting the theory of the Orlicz spaces. 
By virtue of these a priori estimates we then argue in a manner similar to that in [Arch. Rational 
Mech. Anal. 173 (2004), 297-343] to pass to the limit and show that the limiting functions are 
indeed a weak solution which satisfies the mass and momentum equations in the entire space-time 
I domain in the sense of distributions, and the energy equation in any compact subset of the "fluid 
! region" . 

^ ' Keywords: Global weak solutions, 2D Navier-Stokes equations, heat-conducting flows, 

■ spherically symmetric solutions, Orlicz spaces. 
^ 

'nI" 1. Introduction 

o ■ 

^ . The two-dimensional Navier-Stokes equations of compressible heat-conducting fluids express 
the conservation of mass, and the balance of momentum and energy, which can be written as 
' follows in Eulerian coordinates. 

ft + div(H = o, (1.1) 

(^u)i + div(^u ® u) + yP{Q, 6) = fiAu + (A + /x) Vdivu + ^f, (1.2) 
{gE)t + dw{gE + P{g, 9)u) = a(^k9 + ^ii\u\^^ + fidiv [(Vu)u] 

+Adiv[(divu)u] + ^u- f. (1.3) 

Here g, u = (mi,M2) € K^, and 9, 



\u 
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E = ^ + 9 and P = Kg9 

are the density, velocity, temperature, total energy density and pressure of an ideal gas (with 
unit specific heat), respectively; /i and A are the constant viscosity coefficients satisfying /i > 
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and // + A>0, K>Ois the heat-conduction coefficient; f = (/i, /2) is the external force, and Vu 
denotes the gradient of the velocity vector with respect to the spatial variable x e M^. 

We consider an initial boundary value problem for the system fll.ll) - fll.3p in a ball f2 := {x G 
M^; |x| < i?} with boundary conditions 

d6 

u = 0, — = OondQ (1.4) 
an 

and initial conditions 

ig,u,9)\t=o = {g,uo,9o), (1.5) 

where n is the outer normal vector to dQ. 
In the spherically symmetric case, namely, 

g{t,^) = git,r), u(t, x) = ^(t, r)^ ^(t, x) = 0(t, r), f (t, x) = /(t, r)^, (1.6) 
where r = r(x) := |x|, the system fll.ip - fll.3p takes the form: 

gt + {gu)r + — = 0, (1.7) 
r 

{gu)t + {gu%. + — + Prig, e)-v (ur + -) = Qf, (1-8) 
r \ r / r 

ige)t + igue)r + (e„ + ^ V ^(f?- ^) i^r + -) = Q, (1.9) 

r \ ^ J V r/ 



where 



i^:=A + 2/i, Q = I, i^ur + -j -^dru'^>0. (1.10) 
The boundary and initial conditions become 

u = er = atr = R, (1.11) 

{g, u, 6) \t=o = {q, Uo, 9o), < r < R. (1.12) 

The main purpose of this paper is to prove the global existence of weak solutions to the 
problem f ll.ip -f lT3]) when the initial data and external forces are large and spherically symmetric. 
Our work is motivated by the paper of Hoff and Jenssen ^ where they studied the spherically 
and cylindrically symmetric nonbarotropic flows with large data and forces, and established the 
global existence of weak solutions to the compressible nonbarotropic Navier- Stokes equations 
in the "fluid region". In the entire space-time domain, however, the momentum equation in 
jol, Thorem 1.1] only holds weakly with a nonstandard interpretation of the viscosity terms as 
distributions. A natural question is to ask whether the momentum equation holds in the standard 



sense of distributions. A positive answer was given recently by Zhang, Jiang and Xie [17[ for 
a screw pinch model arisen from plasma physics when the heat conductivity n satisfles certain 
growth conditions. In the present paper, based on some new uniform global estimates of u and 6 
(see Lemma 12.41) which are established by applying the theory of the Orlicz spaces, we can give a 
positive result for the two-dimensional Navier-Stokes equations of compressible heat-conducting 
fluids f ll.ip -f fOj) . improving therefore the result of joj. 

We now give the precise statement of our assumptions and results. The external force f is 
assumed to satisfy 

f G L\{0,T),L°^{n))nL°°{iO,T),L^{n)) (1.13) 
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for each T > 0. The initial data (^o, Uq, 6*0) are assumed to satisfy 

<Qo< Co, < Oo a.e. in fi, (1.14) 

f?0'S'(f?o, uo, 6'o)dx < Co for some positive constant Co, (1-15) 
where S is the entropy density in the form of 

Sig, u, 6) = K<f{g-') + ^{6) + ^\u\' (1.16) 

with \l/(s) = s — logs — 1. We point out here that there are no smallness or regularity conditions 
imposed on f and {go, uo, 6*0). 

Under the conditions f ll.13p -f ll.16p . we shall prove the following existence theorem on spher- 
ically symmetric solutions to the problem f ll.ip -f fT3|l . 

Theorem 1.1. Assume that the initial data (gQ,UQ,6Q) and the external force f are spherically 
symmetric and satisfy the conditions U.l^) - D.15\) . Then, the initial boundary problem 
U.5\) has a global weak solution (g,u,6) in the form of U.6\) satisfying the following: 

(a) The support of g is bounded on the left by a Holder curve r{t) G Cj^'^^ : [0, 00) — > [0, 00). 
Moreover, if J-" is the "fluid region", defined by 

T := {{t,x.) I t > and r{t) < r{t) < R}, 

then n {t > 0} n {r < R} is open set. 

(b) The density g G LJ^^(J^), u and 9 are locally Holder continuous in H {t > 0}, and the 
Navier-Stokes equations / fr7]) -/ fr^) hold in V'{F n {t > 0} fl {r < R}). 

(c) The density g G C([0, 00), W^^'°°(i7)*). Also, g(t,-) = in Q\^, and if gu is taken to 
be zero in VL\^ , then the weak form of the mass equation holds for test functions 
^ e C^{\txM X H).- 



/ ^Vdx 
Jn 



*2 fti 

tl 



[ [ {g-^t + gu ■ VV')dxdt. (1.17) 
Jt, Jq 



(d) The velocity u G Lf^^{[0,oo),W^^^/%n)). For ti < t2, if ip G C\[ti,t2] x ^]) vanishes on 
dQ, then, for i = 1,2, 



guitpd'x 



n 



t2 1-12 

{guitpt + guiU ■ Vip + P{g, 0)4jx,) dxdt 



rt2 p /•t2 p 

/ gfidxdt- / [(A + /i)u^, ■ W + /iVMiV^Aldxdt (1.18) 
Jti Jn Jti Jn 

(e) The gradient V6eL\^^{J^) , and the weak form of the energy equation M.3\) holds for test 
functions ip G C^([ti, ^2] x ^) for which there is anrj > such that supp'ip{t, ■) C {x | r{t) + 
< r(x)} for each t G [ti,t2]-' 

p t2 pt2 P I't2 r 

I gEijjdx - / {gEijjt + {gE + P{g, e))u ■ V^)dxdt = / gi ■ mpd^dt 
Jn u Jti Jn Jti Jn 



t2 



tl Jn 



tl Jn Jti Jn 

kVO + ^/iV|up + /i(divu)u + A(Vu)u j ■ VV'dxdt. (1.19) 
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(f ) The total energy energy, minus the mechanical work done by the external force, is weakly 
nonincreasing in time. That is, if 



e{t,^) + ^\n\' 



dxdt, 



Sit) := / f?(t,x) 
Jn 

then ^ 

S{t) = S{0) + / / ■ fdxdt - lim lim / {^E^){t, x)dx 

Jo Jn ^^^1-^°° Jej<r{^)<b 

as a function of t in P'(0, oo), where = |u-' p/2 + 6^ . 

Remark 1.1. In Hoff and Jenssen proved that ([L2]) holds in V'{J'n{t > 0} n {r < R}). 
Here, by exploiting the theory of the Orlicz spaces we are able to derive some new uniform global 
integrability of the approximate solutions (cf. Lemma \2A\\ to show that (11.21) holds in the entire 
space-time domain in the weak sense (i.e., (d) of Theorem 1.1). Furthermore, this method can be 
applied to the screw pinch model with positive constant heat-conduction coefficient in [l7| and 
the cylindrically symmetric rotating model of (ll.ip -f lL3|) (that is, in the symmetric equations 
(6)-(10) in [ol, we take u = v, w = 0, fi = f2 and /s = 0) to obtain similar results. 

Combining the global a priori estimates derived in Subsection 2.3, we shall prove Theorem 
II. II in Section 4 by the convergence argument similar to that in [9^. For this purpose, we consider 
the approximate solutions {g^,u\6^) of the problem (ll.ip - (ll.5l) in the annular regions : = 
{x I Ej < r(x) < R}, where Ej is a sequence of positive inner radii tending to 0. Since the 1/r 
singularity in the equations f ll.7p - (ll.9p plays no role at the stage when Ej is fixed and positive, 
the global existence of approximate solutions {g^,u^,9^) for fll.ip - (ll.5p can thus be shown in a 



manner similar to that in [9|, Il7| . However, to pass to the limit as j — )■ 00 and to show the global 
existence of weak solutions to the original problem (ll.ip - (ll.5p . we need some £j -independent a 
priori estimates. This will be done in Sections 2 and 3. We first prove the global estimates in 
Section 2, where we derive the standard energy-entropy estimates in Subsection 2.1, and apply 
these estimates to establish a new uniform integrability of the approximate solutions in the entire 
spacetime domain by exploiting the theory of Orlicz spaces in Subsections 2.2 and 2.3, which 
is crucial in the proof of (d) of Theorem 1.1. Then, in Section 3, we list the well-known the 
pointwise bounds for qj and 9j as consequences of the energy and entropy estimates. These 
pointwise bounds are independent of Ej, but only away from the origin of Lagrangian space. 



More precisely, as in [9|], for any given h > we define the particle position rj^{t) by 

(f{t, r)rdr 

from which and the standard energy-entropy estimates it follows that there exists a positive 
constant C{h), depending only on h > 0, such that r;^(t) > C{h) > 0. With this observation, we 
can obtain that for any fixed h > and T > 0, there is a positive constant C(T, h), depending 
only on h, T and the initial data, such that 

C(T, h)-^ < g\t, x) < C(T, h) for any (t, x) G [0, T] x [r^(t), R]. 

Applying these pointwise bounds, we can get a number of higher-order energy estimates for the 
approximate solutions in Subsection 3.2, which are also independent of Ej and only away from the 
origin of Lagrangian space. These £j-independent bounds enable us to define the "fiuid region" 
J-" (see (a) of Theorem 1.1) and to obtain the uniform Holder continuity of the quantities on the 
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compact subsets of J-" fl {t > 0} (see (b) of Theorem 1.1). Finally, all the assertions of (a)-(f) 
indicated in Theorem 1.1 will be proved in Section 4 by the convergence arguments adapted from 
Hoff and Jenssen's paper j^. We note that the final step of this argument provides a sort of a 
posteriori validation that the equations fll.7l) - fll.9p are indeed the correct forms of the general 
system fll.ll) - fll.3p in the symmetric case considered here. 

As pointed out in we still do not have sufficient information to infer that r(t) = 0, nor do 
we know whether solutions exist for which r 7^ 0. The analysis simply shows that r{t) may be 
positive, and that, if it is, a vacuum state of radius r{t) centered at the origin. In any case, the 
total mass is conserved in the spherical case, as is clear from (c) of Theorem 1.1, and the total 
momentum is zero because of symmetry. 

We show in (e) only that the energy equation holds on the support of g, rather than in the 
entire space-time domain (0, 00) x Q. This is partly due to that we cannot obtain higher global 
regularity of 6 and u. We may regard the restriction in (5) that the test function be supported 
in J-' as reasonable, since there is no fluid outside J-", and the model is not really valid there. 
Additionally, the failure of the analysis to detect whether or not energy is lost ((f) of Theorem 
II. ip calls into question the adequacy of the mass, energy, and entropy bounds in Lemma 2.1, 
which are the only known (global) a priori bounds in the multidimensional case now. 

We end this section by mentioning some related existence results for large data in the mul- 



tidimensional case. The global existence of weak solutions was ffist shown by Lions [16j for 
isentropic flows under the assumption that the specific heat ratio 7 > 3n/{n + 2) where n = 2,3 
denotes the spatial dimension. Then, by using the curl-div lemma to delicately derive certain 



compactness, and applying Lions' idea and a technique from [12|, Feireisl, Novotny and Petzel- 



tova 0, Isj extended Lions' existence result to the case 7 > n/2. For any l<7<n/2, a global 
weak solution still exists when the initial data have certain symmetry (e.g., spherical, or axisym- 



metric symmetry), see [8|, [ll|-[l2|. For non-isentropic flows, the global existence for general 
data is still not available. Recently, under certain growth conditions upon the pressure, viscosity 
and heat-conductivity (i.e., radiative gases), Feireisl, et al. obtained the global existence of the 
so-called "variational solutions" in the sense that the energy equation is replaced by an energy 
inequality, see jsf for example. However, this result excludes the case of ideal gases unfortunately. 
The global existence of a solution for large data in the non-isentropic case needs further study. 



2. Global Estimates 

In this section we derive a priori global estimates for any smooth (approximate) solution [g"^, 
u'^, O'^) of f ll.7p -f lTT^ together with additional boundary conditions: 

= Qe = Q at r = e and R. (2.1) 

We assume that the initial data and force are smooth and satisfy the bounds fll.13p - fll.15p 
with constants independent of e and 

/R i-R 
^gdx = Mq:= j ^odx. 

We refer to Section 4 for a brief discussion on the existence of such approximate solutions. As 
discussed in Section 1, we shall eventually take a sequence of inner radii Sj — ?■ to prove Theorem 
1.1. Since £ > is fixed for the time being, we suppress the dependence on j. 

2.1. Energy and Entropy Estimates 

We start with the following lemma which states the standard energy and entropy estimates 
for these approximate solutions. 
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Lemma 2.1. Let ((f , , 9^) he a smooth solution of l[1.1/\j - (rj^) defined on [0,T] x [e, -R] with 
boundary conditions 1^2. Then, there are constants Mq and C(T), such that 



R 



g^{t, r)rdr = Mq, 



{g'E'){t,r)rdr <C{T), 



[g'S'){t,r)rdr + 



t rR 



Je 



(2.2) 
(2.3) 

rdr < C(T) for all t G [0, T], (2.4) 



where E'^ = {u'^Y /2 + 9^ , is the entropy density defined in M.10\] and Q is given in U.10\) . 



Proof. The bounds f l2.2p -f l27i|) are the standard energy estimates which follow directly from 
the equations f ll.7p -f lL9|) . the boundary conditions and the assumption f ll.lSp on the external 
force. □ 

2.2. Excursion to Theory of the Orlicz Spaces 

Before deriving the global estimates on the temperature and velocity, we recall some well- 
known results concerning the Orlicz spaces (see, for example, [H, [l5[ for details), which are often 
used to investigate the 2D compressible Navier-Stokes equations (see 13|, llO|, l2| for example). 



Definition 2.1 (Young's function). We say that ^ is a Young's (or N-) function if 

$(t) = I </)(s)ds, t > 0, 
Jo 

where the real-valued function defined on [0, oo) has the following properties 

0(0) = 0, 0(s) > 0, s > 0, lim^^oo 0(s) = oo, 
is right continuous and nodecreasing on [0,oo). 



We define 



i){t) = sup s, t>0, ^(t) 

{<P(s)<t} 



Then is a Young's function as well. We call \E' the complementary Young's function to $. // 
$ is complimentary to , then \l/ is complimentary to 

Definition 2.2 (Orlicz spaces). Let fl be a domain in M" and let ^ be a Young function. The 
Orlicz class Kc^iVL) is the set of all (equivalent classes modulo equality a.e. in VL of) measure 
functions u defined on f2 that satisfy /^^ $(|M(x)|)dx < oo. The Orlicz space L$(f2) is the linear 
hull of the Orlicz class Kq,[VL), that is, the smallest vector space that contains Kq,[VL). The 
functional 



inf jfc > 



MX 



k 



dx < 1 ^ < OO 



is a norm on L$(f2). It is called the Luxembourg norm. Thus, L<^{VL) is a Banach space with 
respect to the Luxembourg norm. 
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Definition 2.3 (Cone condition). Let y be a nonzero vector in R"'. Let Z(x, y) be the angle 
between the position vector x and y. For given such y,h>0, and k satisfying < k < it , the 
set 

A = {x e I X = or < |x| < /i, Z(x, y) < k/2} 

is called a finite cone of height h, axis direction y and aperture angle k with vertex at the origin. 

Q C M" satisfies the cone condition if there exists a finite cone A, such that each x G zs the 
vertex of a finite cone Ax contained in Q and congruent to A. 



Now, we define 

M = M{s) := (1 + s)ln(l + s) - s, N = N{s) := e' 



1, H = His) :-- 



1. 



Next, we list some basic facts on tlie Orlicz spaces Lm{^), Lj\f{Q) and Lh{^)- 

(a) M and N are tlie complementary Young's functions (see [l|, 8.3]). 

(b) Let il(x.) G Lm{^) and f (x) G Lm{^)- By virtue of the generalized Holder's inequality 
0, 8.11]), we have uv G L^{il) and 



see 



uvd'K 



< 2\\u\ 



M{n)\\v\\N{n) 



(2.5) 



(c) Let f2 be bounded and satisfy the cone condition in R^. By virtue of [H, Theorem 8.12, 
8.25 and 8.27], we have for any p > 1 that 



W^-\n) ^ Lh{^) ^ Ln{^) ^ LP{n) and w^'\n) LN{n). 



(2.6) 



(d) Denote by En{Q) the closure of the set of all bounded measurable functions on Q with 
respect to the Luxembourg norm || ■ \\N(n)- Then, the Orlicz space Lm{^) is -Eiv-weakly 
compact, i.e., for any sequence {f„} G Lm{^) uniformly bounded, there is a subsequence 
of {vn}, still denoted by and a t> G Lm{^), such that 



/ fnV^dx — / vifdx for any ip G E]\i{il) 
Jn Jn 



see 



131, 3. Appendix]). 



2.3. Global Estimates on the Temperature and Velocity 

Now, we are in a position derive the global estimates on temperature and velocity. First, we 
define 



^~^(-,|x|) 



^^(■, |x|), £ < |x| < R, 
0, < |x| < £, 



^~^(-,|x|) 



r(-, |x|), £ < |x| < R, 

r(-,£), < |x| < £, 



and make use of (I2.2l) - fl2.4l) and the definition of Luxemburg norm || ■ ]|L^f(n) to deduce that there 
exists a constant Ci(T), such that 



0^{t, \x\)\\mn) = Mo, Wg^it, \x\)\\LM{n) < Ci{T), Vt G (0,T), 
^~-(t,|x|)lnfl + ^~'^(t,|x|))dx<Ci(T), Vt G (0,T), 



and 



Jn 



Vln ( 1 + (t, x) ) dxdt < Ci (T) . 
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(2.7) 
(2.8) 

(2.9) 



Notice that 6^{t,r) is a smooth function in {e,R), we can easily verify that 

ln(l + e^{t, |x|)) G L\0, T; W^'^Q)). (2.10) 

Furthermore, by (12. TP and Lm L^{^)i we infer that there exists a constant C2(f2), such that 

C2{n) < \\g-{t, \^\)\\L,m for any t G (0,T). (2.11) 

At this stage we shall need two auxiliary results: 1) The first one is a revised version gener- 
alized Korn-Poincare inequality (see jsi Theorem 10.17]) in the case of the Orilcz spaces. The 
idea of the proof is essentially the same as that used in the proof of js], Lemma 3.2] under triv- 
ial modification. 2) The other one is the revised version Sobolev embedding in two dimensions 
which will be used to derive bounds of the temperature. This is an idea due to Lions who ever 
used similar embedding to derive the global integrability of the temperature in the proof of the 
existence of weak solutions to the stationary problems for the full compressible Navier-Stokes 
equations in a bounded domain C 16|, (6.204) in Section 6.11]. These two auxiliary results 
are formulated in the following two lemmas. 

Lemma 2.2 (Generalized Korn-Poincare inequality). Let Q G M."^ be a bounded domain satisfying 
the cone condition. Assume that v G W^''^{Q), and g > satisfies 

o<Ci< \\g\\L,m<C2. (2.12) 

Then there is a constant C3 depending solely on Ci and C2, such that 

||t;|L.(^)<C3(Ci,C2) 



Proof. We prove the lemma by contradiction. Suppose that the conclusion of Lemma 2.2 be 
false, then there would be a sequence {^?n}^^Ll of non-negative functions satisfying fl2.12p and a 
sequence {vn}'^=i C W^''^{il), such that 



(^||Vi;„||l2(q) + J ^„|t;„|dx^ , C„ ^ +00. (2.13) 



Setting Wn = Vn\\vn\\ j^^^y making use of (12. 6 p and f l2.13p . we find that 



Wn ^ w = —1= strongly in Lf^i^C). (2-14) 



In view of the hypothesis fl2.12p . we see that 

/ f?nV^dx — )• / ^(/)dx for any Lp G En{^)- (2.15) 
Jn Jn 

Thus, by virtue of (IZMl) . ( l2J[5l) and (ES]), one has 

lim / {QnWn — Qw) dx = lim / Qn{wn — w)dx + lim / {Qn — Q)w(h^ = 0. (2.16) 
The identity fl^TTB]) . together with flZT^ and (EUD, yields 

lim / QnWndx. = / gwdx. > 0. (2-17) 
On the other hand, fl2.13p implies 

lim / gnWnd:>c = 0, (2-18) 

which contradicts with fl2.17p . Therefore, the conclusion of Lemma 2.2 remains true. □ 
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Lemma 2.3 (Sobolev embedding). Let Q be a bounded domain in and 

Q = {e I ||ln(l + 0)||^,.(^)<Ci}. (2.19) 

Then, for any q > 1, there is a constant C2 depending solely on q, Ci and Vl, such that 

||^IU.{f^)<C2(g,Ci,n) for alio eQ. 

Proof. We use (ICTj) and (|Ml) to infer that 

A := ||ln (1 + e)ll^^(^) < C,{n) ||ln (1 + ^)|UM(n) < C,C,. 

Easily, it suffices to consider the case A 7^ 0. By the definition of the Luxemburg norm || ■ ||Lj^(r2), 
we obtain 

j jexp [A-iln(l + e)]^-ljdx< 1, 

1 (l + 0)A''i>-(i+^)dx<l + |(]|. 
Jo. 

[ e''dx< [ ^Mx+ [ (l + 0)Mx 

<(e''^' - iy\Q\ + [ (1 + e)^"'"^(^+'')dx 



which yields 
Hence, 



< 



This completes the proof of Lemma 12.31 □ 
Thus, with the help of Lemma [2.31 and the estimates fl2.7p - fl2.lip . we conclude 
||ln(l + e-it, \M)\\LH(o,nw^^Hn)) < C^{T,Mo,n), 
which together with Lemma 12.31 gives 

Wit, \x\)\\LHio,nL.m) < Wit, |x|)|U2((o,T),L.(n)) < CiiT,Mo,q,n) (2.20) 

for any q G [1,cxd), where fi^ = {x G | £ < |x| < R}. Furthermore, using f l2.4p . fl2.2Up and 
Holder's inequality, we get 

^ (2-21) 

rdrdt) Cy\T,Mo,2,Q) ^ ' 



< 




lo Je i0^nt,r) J 
<C,iT,Mo,n). 

Now, we denote u^(t, x) := u'^it, |x|)x/|x|, which is a smooth function in Q"^. By a simple 
calculation, we find that 

iVu^P = «)2 + (u^)C + < 2^/2|Vu^|^ for any C > 1. 



If we utilize (12. 4p . f l2.20p with g = 2 and Holder's inequality, we obtain 



]£\2 



(t, r)rdrdt 



1/3 



UO Je 

< Ce{T,Mo,n). (2.22) 

Since u^(-,x)|gj7e = 0, we extend by zero outside and employ Sobolev's inequality and 
fl2:22|) to deduce that 

^ (^j |u^|Mx^ dt<C7{T,Mo,n). (2.23) 

Finally, combining (ESQ]) -(ESS]) with (12:231) . we conclude 

Lemma 2.4 (Global estimates of u and 6*). Under the assumption of Lemma l2J[ there are 
constants Ci and C2, such that 

J {eyit,r)rdr\ dt < Ci(T, Mq, g, fi), (2.24) 



T / pR 



T pR 



Je 







^ |ft^(t,r)|rdr^dt <C2(T,Mo,fi), (2.25) 

rdrdt <C2{T,Mo,n), (2.26) 

(M")VdrJ dt<C2iT,Mo,n). (2.27) 



T / /.iJ X 1/3 



3. Local Estimates 



In order to taking to the limit as e — > 0, we will need further uniform bounds of higher 
order derivatives. Such bounds will be obtained away from the origin of Lagrangian space in the 
following sense. Define a curve rf^{t) foTh>0 by 

h= / g%t,r)rdr. (3.1) 



Easily, by ([LTD, 

Thus r|(t) is the position at time t of a fixed fluid particle. Furthermore, an easy estimate, 
based on Jensen's inequality and boundedness of / Q'^'^{Q~^)rdr in ( 12. 4p (see ( 11.16^ ). shows 
that /i — i- at a uniform rate as rf^{t) — )■ 0. That is, given h > 0, there is a positive constant 
C = C{h) independently of e and T, such that 

rlit) > C{h)-\ (3.2) 

Using (13.21) . we can derive pointwise bounds for the approximate density and temperature, 
which are valid away from the origin h = of Lagrangian space, but independent of e. The idea 
of deriving the pointwise boundedness was used first by Kazhikov and Shelukhin 1J| , and later 
adapted by Frid and Shelukhin [3], and by Hoff and Jenssen jij in a nontrivial way to show a 
pointwise boundedness similar to that given in Lemma 13.11 below. 
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Lemma 3.1 (Pointwise bounds). Given h > and T > 0, there is a constant C = {T,h), 
independent of e, such that, if rl{t) is given by 1^3. then 

C-^ < g'it.r) <C for re [rl{t),R] and t e [0,T], 

and ^ 

^ |r(r,-)L,oodr<C 



where || ■ \\h,oo denotes the L°°-norm over [r^ (t) , i?] . 

Proof. Taking n = 2 and m = 1 in the proof of [ol, Lemma 2], we immediately obtain Lemma 
EH □ 

Next, we shall make use of a cut-off function which is convected with the flow and vanishes 
near the origin. The cut-off function is constructed as follows: For given e and h, we can fix a 
smooth, increasing function 4>o{r) with (po{r) = on [0,r/j(0)], < 4>o{r) < 1 on (rft(O), 2r/i(0)) 
and (po{^) = 1 on [2rh{0),R], and then define 4>(t,r) to be the solution of the problem 

(Pt + u^r = 0, 0(O,r) = 0o(r). (3.3) 

We choose (po so that 

<Po{r) < C(/i)(/)[f"^^/^ for some p > 2. 
Thus, we can easily show that this boundedness persists for all time, i.e., 

\Mt,rm < CiT,h)(l){t,r{t)YP-'yP. (3.4) 

We shall take p so large that the exponent on the right-hand side of (13. 4p is close to one. Notice 
that here we have suppressed the dependence of on e and h. 

As in [9i], we now introduce three functionals of higher-order derivatives for {u^, 9^): 



A{T):= sup a{t) [ (j){t , r) ( u'^ + {t,r)rdr 

0<t<T Jruit) \ ^ J 



+ / a(t)0(t,r)(u^)^(t,r)rdrdt, 

Jo Jrn{t) 

B{T):= sup a{t) / 0(t, r)(r)2(t, r)rdr + / / a{t)^{t,r){e'^f{t,r)rdrdt, 

0<t<T Jr^(t) Jo Jr^t) 

R pT pR 

e\2u , / / „2/.N ,2/, ^\lne\2i 



V{T):= sup a\t) / 0"(t, r)(e^)"(t, r)rdr + / / a'{t)(f)\t,r){ey{t,r)rdrdt, 

0<t<T Jr^it) Jo JrUt) 

where a{t) = min{l,t}, "dot" denotes the convective derivative dt + udr, and we have again 
suppressed the dependence on e and h for simplicity. Thus, we have the following estimates. 

Lemma 3.2 (Higher order boundedness). Let h > and T > be given. Then there is a 
constant C = C{T, h), such that 



[ [ + {t,r)rdrdt<C{T,h). 

Jo Jr^t^t) V ^ / 



and 

AiT), B{T), V{T) < C{T,h). 
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Proof. This lemma can be shown following the same procedure as in the proof of jol, Lemmas 
4 and 6] with taking m = 1, v = and w = 0. We should point out here that in the proof one 
should make use of Lemma 2.1, Lemma 3.1, (13. 3 p and (13. 4p . □ 

As the end of this section, we give some uniform integrability estimates. To describe these, 
we define the strictly increasing, convex function G by 

G : [1, oo) ^ [0, oo), G{y) := y\ogy. 

Then G^^ : [0, oo) [1, oo), and one can define for r, c > the function 

u{r,c) ■.= r + rG-^ (^^y (3.5) 

It is easy to see that for each fixed c the function r uj{r,c) is continuous and increasing on 
(0, oo), and that 

lima;(r, c) = 0. 

r-s>0 

Finally, ii E G [0,i?], we define l-E] := J^rdr. 

Lemma 3.3 (Uniform integrability). Let u be the same as in Ii3.5\) . 

(a) Given b > and T > 0, there is a constant G = G{T, b), such that 





u 







b,oo 



+ ||log(^^Vl)||,,oo dt<C(T,6) foranytiMe [0,T], 



(3.6) 



where 9\/ 1 = max{0, 1}. 
(b) If e > 0, and g : [e, i?] — t- M is strictly positive and satisfies 



R 



Q \ogQ rdr 



< G. 



Then, for any measure set E C [e, R\, 

g rdr 



<coi\E\,G). 



(c) Let b > and T > 0. Then there is a constant G = G{T, b) such that, if fort G [0, T], E{t) 
is a measurable subset of [&, -R], and if {g,u,9) = (f)"^, u^, 6*^) is the approximate solution 
described at the beginning of Section 2, then 



Irdrdt < oo 



JE{t) 



grdrdt, G{T, b) 

[ [ gu^rdrdt < G{T, b) u ( [ [ grdrdt, G{T. 

Jo J E(i) \Jo JeH) 



T 

J E{t) 
T 



1/4 



Proof. The estimate ( 13. 6p is a consequence of the entropy estimate (12. 4p . The uniform inte- 
grability bounds in (b) and (c), which are important in showing the limit of these approximate 
solutions to be a weak solution in Section 4, can be obtained using arguments similar to those 
used in j^. Lemmas 8 and 9] , and hence we omit the proof here. □ 
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4. Proof of Theorem 11.11 



By virtue of the a priori estimates derived in Sections 2 and 3, we are now able to prove our 
main theorem by taking appropriate hmits in a manner analogous to that in 

To begin with, we denote by a standard mollifier (in r) of width 5, and for e > 5 we define 
the smooth approximate initial data {qq^ ,Uq^ ,9q^) to {qo,Uq,9q) as follows: 

1) Extend qq by its average value outside [e,-R], mollify with if^, restrict to [e,/?], and then 
multiply by a constant to normalize the total mass to be Mq = Qordr. The resulting density 
function is denoted by QQ^{r). 

2) Redefine Uq to be zero on [0,2e] and [R - 25, R], then mollify with Hs to get Mq . Note 
that Uq^ is identically zero in a neighborhood oi r = e oi R. 

3) Redefine 9q to be its average value on [0, 2e\ and [R — 25], then mollify with Hs to get 6q^. 
Note that is constant in a neighborhood of r = £ or i?. 

The resulting data , Uq^, 9q^ then satisfy the hypotheses fll.l4p and (11.151) with the con- 
stants independent of e and 6. Thus, there is a global-in-time smooth solution (f)^'*^, u^'*^, 9'^'^) of 
the system (ll.7p - (ll.9l) with the initial boundary conditions (II. 121) and (12. ip . This is a result of 
Frid and Shelukhin's work [7] in the annular domain {e,R). Next, we want to pass to the limit 
to get a global weak solution. As in (13. ip . we define the particle path r^''^(t) associated with the 
approximate solution {g'^'^,u'^'^,9'^'^) by 

/•'■;/(*) , 

h= g''\t,r)rdr, h,e,6>0. (4.1) 



4.I. Convergence of the Approximate Solutions 

By the a priori estimates established in Lemmas I3.1H3.2[ 12.11 and 12.41 we have the following 
three propositions, which imply that there is a subsequence {£j,Sj) — )■ (0,0), such that the 
approximate solutions and their associated particle paths are convergent. 



Proposition 4.1. Let (g^'^,u'^'^,9^'^) andr^^^{t) he as described above 



(a) There is a subsequence (ej,5j) — )■ (0,0), such thatr'ff{t) converges uniformly for (t,h) 
any compact subset of [0, oo) x (0, Mq), and the limit rh{t) is Holder- continuous in (t, h) 
on any compact set. 

(b) Ifr{t) ■= \imh^orh{t), then r{t) G C'^^l^^l^^ ^'^^ ^im^odt) = 0. 

(c) // the "fluid region" T is defined by 

J= ■= {(t,r) I r{t) <r <R, 0<t< oo}, 
then J-" n {t > 0} n {r < R} is an open set. 



Proposition 4.2. Let the hypotheses of Proposition 4-i be satisfied. Then there is a further 



subsequence, still denoted by {ej,6j), and limiting functions u and 9, such that 

uniformly on any compact subset of (1 {t > 0}. The functions u and 9 are Holder- continuous 
on any compact set. Furthermore, for any T > 0, 



m"^'^^ ^ u weakly m L^/%{0, T), Wi;^^\0, R]) (4.2) 
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Proposition 4.3. Assume that the hypotheses of Proposition 4-2 hold. Then there is a further 
subsequence (ej, 6j) — (0,0) and a function g{t,r) such that 

g'^'^'it,-) g{t,-) m H'^{\r{t) + r],R],rdr) 

and 

g'^'^'it,-) g{t,-) m L^([r(t) + r^, rdr) 
for all t G [0, T] and all rj > 0. In addition, if g'^^'^^{t, ■) is taken to be zero for r < ej, then 

g'"^'it, ■) -> zn L^([0,r(t)],rdr) when r{t) > 0. 

Also, for h > and T > 0, there is constant C = C{T, h), such that 

C'\T, h)< g< C{T, h) forO<t<T and rhit) < r < R. 

Finally, for h > and t >0, 

rrhit) 

h= g{t,r)rdr. (4.3) 

Jr(t) 

By virtue of the a priori estimates established in Sections 2 and 3, one can show Propositions 
4.1-4.3 in the same manner as that in the proof of j^, Propositions 1-3], except (14.21) of Proposition 
4.2 which are obtained by applying the uniform global estimates given in Lemma [231 In addition, 
the identity (14. 3 p shows that mass is conserved for the limiting solution 

i-R rR 

Mq = / g{t,r)rdr = / gQ{r)rdr. (4.4) 

Jr{t) Jo 

4-2. Weak Forms of the Navier- Stokes Equations 

We now turn to the proof that the limiting functions are indeed a weak solution of the 
Navier-Stokes equations in [0, oo) x in the sense of Theorem 1.1. 

First, the limiting functions g, u and 9 have been defined in the fluid region J-" but not 
elsewhere. We therefore define g, gu and gO to be identically zero in the vacuum region T'^. As 
in Section 1 we let r(x) = |x| and define the velocity vector u : [0, oo) x by 

u(t,x) = M(t,r)-. (4.5) 
r 

Abusing notation slightly, we also write f)(t, x) and 6'(t, x) in place of f)(t, r(x)) and ^(t, r(x)). 
Similar notation applies to the approximate solutions, for which we now write in place of u^^'^^ , 
etc. 

We first show that [g, u, 6) satisfies the weak form (I1.17P of the mass equation. 
Proposition 4.4. Let (g, u, 6) be the limit described above in Propositions \4.1\ -4-3. Then, 

(a) The weak form ( [i.i7| ) of the mass equation holds for any test function : [ti, ^2] x ^ — 
R. 

(b) g e C([0,oo),iyi'°°(fi))*. 

(c) g^/^uE L°°{[0,oo),L^in)). 

(d) ue L^/^{0,T;W^'^/^{n)). 
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Proof. The assertions (a)-(c) follow from the analogous arguments as in the proof of jol, Propo- 
sition 4], where we have made use of Lemma 3.3 and Proposition 14. ltH73l 

To show (d), we use (12. 26 p . f l2.27p . (14. 2 p and the lower semi-continuity to deduce that 



T i-R 

Jo 



u. 



1 4/3 



+ 



and 



u 
r 



4/3 



rdrdt < lim lim 

b— s>0 J— 5>oo 



J I^J u^rdrj dt<limhm J (^J w^rdrj dt<C{T). 



T pR 



i|4/3 



+ 



r 

1/3 



4/3N 



rdrdt < C(T) (4.6) 



We can compute that 

a 



w,(t,|x|)^ + «(t,|x|) 



,l<^,J<2, 



(4.7) 



thus, |Vu|4/3 < 2^/3 



14/3 



(Jl)). This completes the proof. 



+ |M/r|^/3) and jup = u\ Hence, we find that u G L^/^{0,T;W^'^/^ 



□ 



Next, we show that the weak form of the momentum equations in the spherically symmetric 
case holds. 

Lemma 4.1. Let g, u and 9 he the functions given in Propositions 4-2 and 4-3. Let ti < t2 and 
(j) be a C^-function on [^1,^2] x [0,-R], such that 4>(t,0) = (p(t,R) = for t G [ti,)f:2]- Then, the 
following identity holds. 



R 

QU(j)rdr 
t2 rR 



t2 



t2 rR 



t-i Jo 



rdrdt — u 



QU(j)t + QU (pr + P{q, I 
t2 rR 



\ + - 
r 



rdrdt 



u 

Mr H — 



(4i 



_l — rdrdt. 



Proof. We first consider a simpler case in which the test function vanishes in a neighborhood 
of the origin. Assume that ip is a. C^-function on [ti, t2] x [0, B\ satisfying = on [0, h] for some 
6 > 0. Then, applying the proof of 0, Lemma 10 (a)] and combining with the weak convergence 
(14.21) . we can easily show that the weak form of the momentum equation (II. 8p holds for the test 
function ip: 



R 

guiprdr 
t2 rR 



t2 



ii -'0 



QUIpt + QU^Ipr + P \ i!r + 



rdrdt 



ti -'0 



gfiprdrdt — v 



t2 rR 



ti JO 



U 



Ur 



(4.9) 



rdrdt. 



To extend the identity (14. 9 p to the case that test functions are supported in [0, R], we fix an 
increasing function x • [0, 00) — ?■ [0, 1] with x = on [0, 1] and x = 1 on [2, 00), and define 
x\r) := xir/b) for b > 0. Let be a function on [ti, ta] x [0, R] such that 0(t, 0) = 0(t, R) = 
for t G [ti,t2], and define 0^ := Then, the previous lemma applies to the test functions 

0'' = x''4'- We obtain 



R 

Qwy^cprdr 

t2 rR 



t2 



tl 



t2 rR 



tl JO 



tl Jo 



gfx (prdrdt — v 



t2 rR 



tl Jo 



U 

Ur + - 

r 



(xV). 



+ 



rdrdt. 



rdrdt 



(4.10) 
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The first, second, fifth and sixth terms in f l4.10p converge to the corresponding terms in (14 .Sp as 
6 — )■ by virtue of the dominated convergence theorem. For the third term we have 

I ^M^(xV)rrdrdt = / / ^^^(x^^ + xVr)rdrdt, (4.11) 
ti Jo Jh Jo 

and the second term on the right-hand side of (l4.1Up clearly tends to the third term in (14. 8p 
as 6 — !■ 0. Since 4>(t,0) = 0, we can write (p(t,r) = rip(t,r) for some smooth ip. Thus, due to 
|<9rX^I < C/6 we can bound the first term on the right-hand side of (14. lip by 

/"<2 rb ^ /•t2 /•2b 

/ / gu^C-\ip{r)\rdrdt <C / gu^drdt, (4.12) 

Jti Jo " Jti Jb 

which goes to zero as 6 — )■ by utihzing the boundedness of the hmiting energy. Moreover, the 
same argument apphes to the fourth term in (l4.1Up . Finally, for the last term on the right-hand 
side of (l4.1Up . we have 

+ -) (^l + -] rdrdt = r r (ur + -) f 0Xr + + — ^ ^^drdt. (4.13) 



r/\ 1" I I I \ iV 

.1 - u ^ V ' / Jti Jo ^ ' 



Similarly to (I4.12p . we deduce that 



^ ' ^ {ur + ^) 0Xr^drdt <C (ur + ^) rdrdt, (4.14) 

which tends to zero as 6 —t- by ( 14. 6p . On the other hand, letting 6 — )► in (14.13p . using the 
dominated convergence theorem and (I4.14p . we conclude 



This completes the proof of the lemma. 

Remark 4.1. Applying (I2.26p -f l2.27p and the proof of Lemma [4. II to jij. Lemma 11], we can see 
that ^ ^ 

hrn^ lim U(j, cj)^) = u + ^ j + rdrdt 

still holds for cylindrically symmetric case in [ol. Lemma 11]. 

Now, we are able to show that the weak form (ll.lSp for the momentum equations in Cartesian 
coordinates is satisfied. 

Proposition 4.5. The weak form ^1.13^) of the momentum equations, as stated in Theorem 1.1 
(d), holds. 

Proof. Given ip as described in the theorem, we define 

(f){t,r) := / ■ip{t,ry)yidSy for fixed i = 1, 2, 



s 
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where 5* C denotes the unit circle. Then, 0) = R) = 0. It thus follows from Lemma 
Othat 

/•H t2 pt2 pU r 

gu(l)t + Qu^(t)T + P + 



gu(f)rdr 
t2 rR 



t2 pR r 

ti Jti Jo 



rdrdt 



ti JO 



gf(f)rdrdt — v 



t2 pR 



ti JO 



Ur + 



U 



(4.15) 



'r H — I rdrdt. 

r 



We convert each of the terms in ( 14.15^ to integrals in Cartesian coordinates involving ip. The 
treatment of the terms involving derivations are very much similar to those in the proof of 
Proposition 4.4 (see [ol, Proposition 4]), except for the last term, which we deal with in details. 
We may rewrite the last term of fl4.15p as 



t2 rR 



ti JO 



Ur + 



U 



t2 rR 



ti JO 
t2 rR 



r -\ — I rdrdt 

- (prr + -[r<p)r 
r J r r 



ti JO 



u 

r ( — 

r 



+ 



t2 rR 
ti Jo 



drdt 

i'xk{t,ry)yiykdSy^ rdrdt 
^a;,(t,ry)dS'y ) rdrdt 



(4.16) 



t2 



ti Jn 

t2 





x|) 




X 





?/'a;^(t, x)rdxdt 



ti Jn 

t2 



M(t, |x|)— V^3.^(t,ry)rdxdt 



ti Jn 



Vui{t, [xDWdxdt, 



where the repeated indexes should be summed. Notice that dx^Uiit, |x|) = dxiUj{t, |x|) (see (14. 7p ). 
we have 



t2 

ti Jn 



Vui{t, |x|) ■ VV'dxdt = / dx^Ui{t, |x|)9x.fcVdxdt 

Jn 



t2 



dx,Uk{t, |x|)9a;^?/'dxdt 



h Jn 

t2 



u Jn 



Mx, ■ VV^dxdt, 



whence, by recalling z/ = A + 2/i, the last term on the right-hand of (I4.15P can be written 

l-t2 r rt2 r 

{X + jj) / Ua,. ■ VV'dxdt + /i / / VwiVV'dxdt. 
Jti Jn Jti Jn 

This completes the proof of Proposition 4.5. 



□ 



Proposition 4.6. The weak form U.19\) of the energy equation, as stated in (e) of Theorem \l.l\ 
holds. The total energy £ of the the limiting functions satisfies (f) of Theorem \l.l[ 

Proof. The proof is the same as that of [9, Propositions 7 and 8], and hence we omit the proof 
here. □ 
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Thus, we have completed the proof of Theorem 1.1. In fact, Part (a) of Theorem 11.11 is just 
Proposition 14 . 1 1 with the semi continuity used imphcitly in the proof. The existence and regularity 
of u, 9 asserted in (b) of Theorem 11.11 follow from Propositions 14.21 and 4.3. The weak forms 
of the mass and momentum equations are proved in Propositions 14.41 and 14. 5[ from which the 
regularity assertions in (c) and (d) of Theorem 11.11 follow immediately. Finally, the results in (e) 
and (f) of the Theorem 11.11 are proved by applying Proposition 4.6. 
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